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ABSTRACT 

In a previous paper we reported a 3(J detection of an absorption line from the Warm-Hot Intergalactic Medium 
(WHIM) using the Chandra and XMM X-ray grating spectra of the blazar H2356-309, the sight-line of which 
intercepts the Sculptor Wall, a large-scale superstructure of galaxies at z 0.03. To verify our initial detection, 
we obtained a deep (500 ks), follow-up exposure of H2356-309 as part of the Cycle-10 Chandra Large Project 
Program. From a joint analysis of the Cycle-10 and previous (Cycle-8) Chandra grating data we detect the 
redshifted O VII WHIM line at a significance level of 3.4cr, a substantial improvement over the l.la level 
reported previously when using only the Cycle-8 data. The significance increases to 4. Oct when the existing 
XMM grating data are included in the analysis, thus confirming at higher significance the existence of the line 
at the redshift of the Sculptor Wall with an equivalent width of 28.5 ± 10.5 mA (90% confidence). We obtain 
a 90% lower limit on the O VII column density of 0.8 x 10'^ cm"^ and a 90% upper limit on the Doppler b 
parameter of 460 km s~'. Assuming the absorber is uniformly distributed throughout the 15 Mpc portion 
of the blazar's sight-line that intercepts the Sculptor Wall, that the O VII column density is w 2 x 10'* cm"^ 
(corresponding to > 150 km "' where the inferred column density is only weakly dependent on b), and that the 
oxygen abundance is 0. 1 solar, we estimate a baryon over-density of ^ 30 for the WHIM, which is consistent 
with the peak of the WHIM mass fraction predicted by cosmological simulations. The clear detection of O VII 
absorption in the Sculptor Wall demonstrates the viability of using current observatories to study WHIM in the 
X-ray absorption spectra of blazars behind known large-scale structures. 

Subject headings: BL Lacertae objects: individual (H 2356-309) — intergalactic medium — quasars: absorp- 
tion lines — large-scale structure of universe — X-rays: diffuse background — X-rays: 
galaxies: clusters 



1. INTRODUCTION 

The total amount of the luminous baryons in the nearby uni- 
verse probed by the stellar light, narrow Lya absorption, as 
well as the X-ray emission from the hot intracluster and intra- 
group medium, accounts for at most 50% of the total baryonic 
matter in the low-redshift universe (e.g., Fukugita et al. 1998). 
Large-scale, cosmological hydrodynamic simulations predict 
that most of the "missing baryons" are distributed as filamen- 
tary structures between galaxies, in the form of a warm-hot 
intergalactic medium (WHIM; T 10^ - 10^ K) (see, e.g., 
Fukugita et al. 1998; Dave et al. 1999; Cen & Ostriker 1999; 
2006; Dave et al. 2001), with typical overdensity 6 between 5 
and 20(Cl 

The first evidence of the WHIM gas came from the de- 
tection of highly ionized, ultraviolet absorption lines in the 
spectra of background quasars. With the Hubble Space Tele- 
scope (HST) and the Far Ultraviolet Spectroscopic Explorer 
(FUSE), a number of intervening O VI absorption systems 
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were detected (see, e.g.. Savage et al. 1998; Shull et al. 1998; 
Tripp & Savage 2000; Oegerle et al. 2000; Tripp et al. 2001; 
Savage et al. 2002; Sembach et al. 2004; Richter et al. 
2004), for the first time revealing the existence of the WHIM 
gas (O VI is sensitive to gas at temperature between lO'* and 
5 X 10^ K). 

While the ultraviolet absorption lines detected in the back- 
ground quasar spectra reveal its existence (see Tripp et 
al. 2006 for a review), the majority of the WHIM is expected 
to be in a temperature region that can only be studied in X- 
rays (e.g., Cen & Osti'iker 1999; Dave et al. 2001). Since the 
launch of Chandra and XMM-Newton nearly a decade ago, 
combined high resolution imaging and spectroscopy has, for 
the first time, made the X-ray study of the WHIM possible. 
The diffuse hot plasma associated with the WHIM can be 
probed in X-ray emission and absorption, with each technique 
possessing certain advantages. The spatial distribution is most 
naturally revealed by imaging the X-ray emission. However, 
because of the low density of the WHIM gas and the contam- 
ination from the Galactic foreground. X-ray emission mea- 
surements are very challenging, and are therefore restricted to 
the highest density tail of the WHIM distribution (e.g., Soltan 
et al. 2002; Kaastra et al. 2003; Finoguenov et al. 2003; Zap- 
pacosta et al. 2005; Lieu & Mittaz 2005; Galeazzi et al. 2007; 
Takei et al. 2007; Werner et al. 2008). For example, the re- 
cent X-ray imaging study by Werner et al. (2008) of a mas- 
sive binary cluster system obtained a 5(t detection of WHIM 
between the clusters with an estimated baryon over-density 
of ~ 150, which is larger than expected for the bulk of the 
WHIM gas (Dave et al. 2001; Cen & Ostriker 2006). 

There are distinct advantages to studying the WHIM in X- 
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Fig. 1 . — Broadband spectrum of obs.#10577 between 5 and 30 A. The red line is a power-law model folded through the detector response; see table 1 for 
model parameters. 



ray absorption (Fang & Canizares 1997; Hellesten et al. 1998; 
Perna & Loeb 1998). In contrast to the X-ray emission which 
depends on p^^^, the column density of WHIM gas revealed 
in the X-ray absorption spectrum of a background quasar de- 
pends on pc;as, which is an advantage for the very low densities 
expected in the WHIM. Moreover, the Galactic foreground 
emission is dwarfed by the continuum of the quasar, while 
any foreground Galactic absorption line is easily separated in 
wavelength from the redshifted WHIM line at cosmological 
distance. However, the spectral resolution and flux sensitiv- 
ity of existing X-ray grating spectrometers severely limit such 
studies. The first study using Chandra did not reveal any de- 
tections because the background quasars were too faint and 
thus did not produce spectra of high enough quality (Fang et 
al. 2001). In fact, most early attempts with Chandra were sim- 
ilarly unsuccessful (e.g.. Fang et al. 2002; Mathur et al. 2003; 
McKernan et al. 2003; Fang et al. 2005). 

The first report of a significant detection of the WHIM in X- 
ray absorption was made by Fang et al. (2002). Using the Low 
Energy Transmission Gratings (LETG) and the Advanced 
CCD Imaging Spectrometer (ACIS) combination in the Chan- 
dra X-ray Telescope, an O VIII absorption line produced by 
the intragroup medium of a small galaxy group along the sight 
line of PKS 2155-304 was detected with more than 4(7 con- 
fidence. This detection later was confirmed with more ob- 
servations using the same instrumental set (Fang et al. 2007), 
along with several other independent analyses (e.g., Williams 
et al. 2005; Yao et al. 2009). However, this line was not de- 
tected when using the Chandra LETG with the high resolu- 
tion camera (HRC), as well as the XMM-Newton X-ray tele- 
scope, although the upper limits from the non-detections are 
fully consistent with the detected line parameters (Williams et 
al. 2005; Fang et al. 2007). Despite the consistency between 
these different measurements, the failure to detect the line 
independently with multiple instruments has generated some 
skepticism (e.g., Richter et al, 2008). Nevertheless, given the 
crucial importance of a firm WHIM detection for our under- 
standing of the cosmic baryon budget and structure formation, 
it is desirable that the result is not contingent only on a single 



instrument (nor, indeed, measurement) which could plausibly 
be affected by hitherto unrecognized systematic effects. 

This point is underlined by the recent controversial claims 
of WHIM detections in X-ray absorption along the sightline 
to the unusually bright blazar, Mkn 421, while in outburst 
(Nicastro et al, 2005). Based on Chandra grating spectra, 
Nicastro et al. reported highly significant detections of two 
WHIM systems (at 3.5(t and 5.8(7), which were not detected 
with a deep XMM-Newton RGS observation (Rasmussen et 
al. 2007). Although the RGS upper limits on the absorption- 
line column densities are consistent with the Chandra mea- 
surements, Kaastra et al. (2006) showed that the line sig- 
nificances reported by Nicastro et al. were inflated because 
insufficient redshift trials were employed to assess the statis- 
tical significance of the absorption systems, given that their 
redshifts were not known a priori. Kaastra et al.'s criticism 
highlights a fundamental limitation with the "blind"-search 
method, which has been the primary tool to search for WHIM 
in absorption. This technique involves observing bright back- 
ground sources at random directions on the sky, in the hope 
that their line of sight intersects a WHIM filament, the redshift 
of which is a priori unknown. 

An alternative strategy is to focus on known foreground 
structures that likely have associated WHIM. The advantage 
of this method is that the absorption signatures can be re- 
garded as highly significant even if they are less prominent 
than those found in a random search, since the redshift of 
the absorber is known a priori from the redshift of the su- 
perstructure traced by its galaxies. Based on cosmological 
simulations of the local universe, Kravtsov et al. (2002) first 
proposed that the WHIM absorption lines are most likely to 
be detected in the local superstructure region. Fujimoto et 
al. (2004) first observed a quasar behind the Virgo Cluster, and 
found low-significance evidence of O VIII absorption from the 
foreground cluster Takei et al. (2007) looked at an active 
galactic nucleus (AGN) behind the Coma Cluster, and also 
marginally detected O VIII and Ne IX absorption from the hot 
gas in the galaxy cluster 

In 2007 we performed such a "targeted" study of the WHIM 
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by observing the blazar, H 2356-309 (z = 0.165), with Chan- 
dra and XMM. The blazar's sight-line passes through multiple 
large-scale galaxy structures, including the Sculptor Wall su- 
perstructure (z ^ 0.03; see Figure 1 of Buote et al. 2009; B09 
hereafter). Blazars are preferred for this type of observation, 
not just because they are X-ray bright, but also because they 
typically lack intrinsic absorption features that may compli- 
cate the search for intervening features. For the first time, we 
detected a WHIM absorption line at the redshift of the Sculp- 
tor Wall with both telescopes, with a joint detection signifi- 
cance of 3cr (B09; specifically, B09 obtained a 99.64% detec- 
tion significance when using the C-statistic and 99.81% when 
using x^)- 

While this result is very promising, in light of the various 
controversies which have plagued past attempts to measure 
the WHIM in absorption, a definitive detection requires fur- 
ther confirmation at a higher statistical significance level. In 
particular, the detection significance was driven mostly by 
the XMM RGS data, with the line only found at 1.1 a in the 
Chandra grating data. To address this, in 2008 we performed 
a deep, follow-up observation of H 2356-309 as part of the 
Chandra Cycle 10 Large Project Program, which we discuss 
in this paper We confirm the presence of the Sculptor WHIM 
line at 3A<t with Chandra alone, while improving the joint 
Chandra-XMM significance to 4-a. This clear detection with 
both satellites constitutes arguably the most compelling X-ray 
evidence to date of a WHlM-like absorber along the line of 
sight to a blazar. 

In this paper we will focus on the Helium-like oxygen Ka 
Une (O VII, rest wavelength 21.6019 A) studied by B09. Oxy- 
gen is the most abundant heavy element in the intergalactic 
medium (IGM), and numerical simulations predict that given 
the WHIM gas temperature, most oxygen should be in the 
form of O VII (Fang et al. 2002; Chen et al. 2003; Cen & Fang 
et al. 2006). Unlike O VI and O VIII, it probes a wide temper- 
ature range of the WHIM gas between 5 x 10^ to 3 x 10^ K. 
Absorption lines from other ion species may also be expected 
but at lower significance, and a companion paper (Buote et 
al. 2010) will discuss those lines. Evidence for the WHIM ab- 
sorption in galaxy superstructures at higher redshifts than the 
Sculptor Wall will be presented in Zappacosta et al. (2009). 

2. OBSERVATION AND DATA EXTRACTION 

H 2356-309 was observed first in 2007 during Chandra Cy- 
cle 8 for 100 ks, and it was observed again in 2008 during 
Cycle 10 in ten separate exposures totaling ^ 500 ks. The 
exposures range from ^ 15 to 100 ks; see Table [T] for de- 
tails. We selected the Chandra Low Energy Transmission 
Grating and the High Resolution Camera (LETG-HRC) com- 
bination. It offers high sensitivity with less instrumental fea- 
tures between 22 and 23 A, where the redshifted O VII line 
is located. The LETG-HRC also offers a constant spectral 
resolution of 0.05 A. The Cycle 10 observation substantially 
improved the count statistics when expressed in terms of the 
counts per resolution element; i.e., the photon counts within a 
bin size of the spectral resolution (50 mA). In the vicinity of 
22 A the combined Chandra data sets from Cycle 8 and Cycle 
10 (excluding sequence 10498, as discussed below) have a 
total CPRE of ^ 460 (including ^100 background counts), 
which greatly exceeds the CPRE of ^ 70 (including ^ 20 
background counts) obtained for just the Cycle 8 data. 

As in B09, we followed the standard procedures to extract 
the spectra. We used the software package ClAO (Version 



TABLE 1 
Chandra OBSERVATION LOG 



Observation ID 


Exposure 
(ks) 


Date 


r 


Flux" 


8120* 


100 


Oct 1 1 2007 


1.896 ±0.033 


1.12 


10577 


82 


Sep 17 2008 


1.810±0.025 


2.08 


10497 


52 


Sep 19 2008 


1.835 ±0.031 


1.95 


10498 


80 


Sep 22 2008 


1.784 ±0.027 


1.94 


10762 


35 


Sep 25 2008 


1.699 ±0.038 


2.04 


10761 


45 


Sep 27 2008 


1.700 ±0.033 


2.20 


10499 


56 


Sep 29 2008 


1.685 ±0.030 


2.01 


10764 


103 


Oct 10 2008 


1.863 ±0.024 


1.74 


10840 


15 


Dec 23 2008 


1.908 ±0.061 


1.87 


10500 


16 


Dec 25 2008 


1.938 ±0.057 


1.89 


10841 


15 


Dec 28 2008 


1.871 ±0.059 


1.95 



b. This sequence was first analyzed in B09. 

4.C0) and calibration database CALDB (Version 3.^ devel- 
oped by the Chandra X-ray Center Most changes in the latest 
release ClAO V4.1 and CALDB V4.1 are related to the Ad- 
vanced CCD Imaging Spectrometer (ACIS) and do not impact 
our analysis. Specifically, we followed the science thread for 
analyzing LETG-HRC dat^3> in which one first generates a 
type 11 pha file, builds the redistribution matrix files (RMFs) 
and the ancillary response files (ARFs), and groups the spec- 
trum. 

As in B09, we have generated our own type II pha file, 
rather than using the file produced by the standard pipeline 
(Reprocessing 111), to take advantage of an improved back- 
ground filter not yet available in the standard processing. This 
LETG background filter recently developed by the Chandra 
X-ray Center (CXC)[]3can allow us to reduce the background 
rate by as much as 50% at the wavelength regions of inter- 
est 20 A) with a negligible loss of source X-ray events 
(^ 1.25%). We therefore reprocessed the data and applied 
this new filter to generate the type 11 pha spectra. 

Unlike ACIS, the HRC does not have the energy resolution 
to sort individual orders, and each spectrum contains contri- 
butions from all the diffraction orders. To account for the 
high-order contributions, it is necessary to build the response 
matrix that includes all the relevant orders. Our experience 
showed that orders higher than the sixth make negligible con- 
tributions. Therefore, as in B09, we ignore them and focus 
only on the first to the sixth orders. For each observation, we 
first built the RMF and ARE for each order; these RMFs and 
ARFs were then combined together, following the procedures 
described on the Chandra websitj^ We also combined the 
negative and positive orders. 

We rebinned each LETG-HRC spectrum to have a min- 
imum of 40 counts per bin. We examined other binning 
schemes (e.g., fixed wavelength bins) and did not find sub- 
stantial improvement in the spectral fitting results. 

Finally, we refer the reader to B09 for details of the XMM 
RGS observation and data preparation. 

3. DATA ANALYSIS 
3.1. Continuum 



^ see http://asc.havard.edu/ciaol 
' see http://asc.havard.edu/caldb' 

see http://asc.havard.edu/ciao/threads/gspec.html 
" see http://cxc.harvard.edu/contrib/letg/GainFilter/software.htmll 



see http://cxc.harvard.edu/cal/Letg/Hrc_QE/ea_index.htniir 
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Fig . 2 . — Chandra data sets, labeled with their observational IDs used in the spectral fitting. The red line in each panel shows the model obtained by simultaneous 
fitting of all the spectra displayed, including the XMMRGS data set. The WHIM absorption line (shadowed region) is located at ~ 22.25 A, and the local Galactic 
line (dashed line) is located at ^ 21.6 A. The shadowed regions show the location of the absorption line, and the size represents the 90% redshift measurement 
error. Note that we experimented with different spectral binning schemes and did not find substantial improvement in the results compared to those obtained 
using our adopted scheme. 



To make an accurate determination of the WHIM absorp- 
tion line properties it is desirable to measure the continuum 
level as close as possible to the line. However, because of the 
contribution from overlapping higher orders diffracted from 
higher energies, as noted above, it is necessary to model the 
continuum over a wide energy band. Consequently, following 
B09 we first fitted a power-law to the broad-band spectrum (5 
- 30 A) for each exposure. (By fitting down to 5 A we account 
for contributions from orders 2-6, including order 3, which is 
the most important higher order.) For determining the WHIM 
(and local absorption) Une parameters, we actually only fit- 
ted the data over a restricted continuum range (2 1-22. 5 A). 
We model the continuum as a powerlaw with its spectral in- 
dex restricted to the 90% confidence range determined from 



the broad-band fit (given in Table 1), but its normalization 
allowed to vary freely. This procedure balances the need to 
account for the contributions from higher orders while allow- 
ing for some adjustment in the local continuum with respect 
to that established by the broadband fit. 

For each observation, we fitted the broadband spectrum (5 
- 30 A) with a power law with the Galactic neutral hydrogen 
absorption. We used the software package XSPECVll.sH 
for spectral fitting throughout the data analysis. Since 
fits of Poissonian data, even when rebinned far more heavily 
than the canonical 20 counts per bin, tend to be significantly 
biased, we performed the fits by minimizing the C-statistic 

See ,http://heasarc.nasa.gov/docs/xanadu/xspec/xspecl l/rndex.html | 
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Fig. 3. — Panel (a) shows the stacked spectra of the ten Chandra data sets, and panel (b) shows the XMM-Newton data set. The Chandra data were stacked 
using an adaptive binning scheme (see text) intended for display puiposes only. The red lines are not produced by fitting the stacked spectrum, but rather are the 
results of simultaneous fitting of the individual Chandra data sets (Figure[2) with the XMM-Newton data set. 



(identical to maximizing the Poisson likelihood function) of 
Cash (1979), which yields less biased best-fitting parameters 
(Humphrey, Liu & Buote 2008)E3. We note that the statistical 
significance does not change appreciably between the and 
C-statistic (see below). 

The column density estimated from the HI map of Dickey 
& Lockman (1990) by the Colden tooQ is Nh = 1.33 x 
10^*^ cm~^, and we fixed the Galactic absorption at this level. 
We also tested our fits with free Nfj and found the fitted Nh 
differed very little, indicating no excess neutral hydrogen ab- 
sorption. Figure [U shows an example of the broadband fitting 
for the obs. #10577, demonstrating that our simple model (red 
line) provides a satisfactory fit to the broadband behavior. 

The photon index and spectral flux of each observation are 
listed in columns 4 and 5 in Table iQ For consistency, we 
also refitted obs. #8 120, listed in the first row, that was ana- 
lyzed in B09, and report its flux including Galactic absorption 



Although it is not strictly necessary to bin the data to obtain the best- 
fitting model when using the C-statistic, it is critical to bin the data in order to 
assess the statistical significance of the line detection (H ump hrey et al, 2009). 
We discuss the impact of different binning schemes in j|4.4| 

see http://asc.harvai'd.edu/toolkit/colden.jsp 

Thi'oughout the paper, all the errors and Hmits are quoted at 90% confi- 
dence level, unless otherwise specified. 



(B09 reported the unabsorbed flux.) The flux is measured 
between 0.5 and 2 keV, in units of 10~" ergs s"'cm~^. The 
flux increased almost by a factor of ^ 2 between the Cycle 
8 (observed in October 2007) and the Cycle 10 (observed in 
late 2008) observations. Overall, the 2008 flux levels are very 
consistent with previous BeppoSAX (Costamante et al. 2001) 
and RXTE (Aharonian et al. 2006) observations. The photon 
index also varied between approximately 1.7 and 2. 

3.2. Absorption Lines 

In Figure|2]we display the spectrum for each individual ex- 
posure over the « 21-22.5 A wavelength range used for our 
analysis. The positions of the O VII Ka absorption lines cor- 
responding to the local (z ^ 0) absorber and the Sculptor Wall 
(z ~ 0.03) previously detected by B09 are indicated. For one 
of the exposures (obs. #10498, not shown), there is also a sig- 
nificant absorption feature located near 22. 05 A, and we defer 
analysis of that data set until discussing systematic errors in 
^4.11 Consequently, the total Cycle 8 and Cycle 10 Chandra 
exposure (omitting #10498) used in our analysis is~ 520 ks. 

The absorption lines are fitted using a model that was devel- 
oped in B09 rather than simply subtracting a Gaussian compo- 
nent from the continuum. This physical, Voigt-broadened line 
absorber model has three free parameters; the column density 
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TABLE 2 
O VII Absorption Line Properties 



Line 


log(N) 
(cm-2) 


b 

(kms-') 


Chandra Redshift 


XMM Redshift 


EW 
(mA) 


Sculptor Wall 


16.84!i;» 


96 « 455) 


0.0306 ±0.0007 


0.0322 ±0.0012 


25.8 ± 10.5 



Local 



17.70!«:54 



-0.0008 ±0.0016 



0.0018 ±0.0018 



17.4± 12.3 



Note. — All the errors are quoted at the 90% confidence level on one interesting parameter (AC = 2.71). The Doppler b parameter was allowed to vary 
between 20—600 km s"' . 



of the absorbing ions, the Doppler-/? parameter, and the red- 
shift. This model is motivated by the absorption line physics, 
and the model parameters reflect directly the physical proper- 
ties of the absorbers. Consequently, this model readily allows 
the line parameters, notably the column density, to be tied 
consistently between the different spectra, unlike ad hoc mod- 
els that simply subtract a gaussian from a continuum model. 
Note following B09, in calculating the optical depth we ig- 
nore the correction term from stimulated emission since it is 
negligible at the gas temperature where the O VII ionization 
fraction peaks. 

To determine the significance of the WHIM line at ^ 22.25 
A, we adopt the following procedure. We first fit the Galactic 
(i.e., local) line in addition to the local continuum (all modi- 
fied by Galactic neutral hydrogen absorption), where, as noted 
above, the latter is restricted to have powerlaw exponent con- 
sistent with the 90% limits that were established by the broad- 
band fit. The redshift of the added line was restricted to match 
the structure in the Sculptor Wall (z=0.028-0.032; B09). We 
then add the WHIM line, fit again, and the resulting C-statistic 
difference between this fit and the previous one (local contin- 
uum -I- the Galactic line) gives the significance of the WHIM 
line. For the Galactic line, we follow the same procedure by 
fitting the WHIM line first and then adding the Galactic line, 
which we restrict to have z<0.0025. 

We have two choices for spectral fitting. One is to stack all 
the observations into one spectrum, and combine the response 
and effective area file of each observation into one response 
file. This approach is advantageous for low signal-to-noise 
spectra. However, due to the variation of each spectrum (both 
power law index and flux), which is typical for blazars, this 
approach will produce significant systematic errors, as de- 
scribed in Rasmussen et al. (2007). So instead of stacking 
all the spectra, we simultaneously fit all the data sets, such 
as we did previously with the Cycle-8 Chandra data and the 
XMM data in B09. Each specti'um has its own RMF and ARF. 
While we allow the continuum parameters (photon index and 
normalization) to vary separately for each observation, the pa- 
rameters for the absorption lines are tied between each Chan- 
dra observation. We also allow the redshift of the XMM data 
to be fitted separately from the Chandra data, which we dis- 
cuss more below. 

In Figure|2]we display for each Chandra spectrum the best- 
fitting model obtained from a simultaneous fit of the Chandra 
and XMM data. Because the absorption lines are not immedi- 
ately apparent to visual inspection in every Chandra spectrum, 
to aid visual interpretation we stacked the individual spec- 
tra into a single combined spectrum using an adaptive bin- 
ning scheme. We show the resulting stacked spectrum for the 
Chandra data in Figure 3 (along with the XMM data), which 
can be considered a lightly smoothed spectrum constructed 




20.0 20.5 21.0 21.5 22.0 22.5 23.0 23.5 
Wnvelength ( A ) 

Fig. 4. — Stacked Chandra data and best-fitting model as in Figure|3]except 
shown for a wider wavelength range. We omit the corresponding expanded 
RGS spectrum because much of the added bandwidth is contaminated by in- 
strumental features. We caution the reader against interpreting uninteresting, 
narrow, statistical fluctuations in the spectra as real features. The likelihood 
of finding such a statistical fluctuation is dramatically enhanced if the fine 
energy is not a priori known, in contrast to the case of the local and Sculptor 
wall features we discuss in this paper 

from some components with a coarser bin scale. Also shown 
are the best-fitting models obtained from the simultaneous fits 
to the individual exposures. We reiterate that the models dis- 
played in Figure 3 are not fitted to the stacked spectrum. In 
Figure |4] we show the expanded Chandra spectrum between 
20.5 and 23.5 A. 

Visual inspection of Figure 3 reveals very clearly the pres- 
ence of the WHIM line near 22.25 A. This is confirmed by the 
improvement in the spectral fits when adding the WHIM line. 
The C-statistic decreases by 14.1 (229 bins) when using only 
the Chandra data, and decreases even more (20.4 for 271 bins) 
when including the XMM data in the fits. The reduction in the 
C-statistic and the total number of bins are over twice that re- 
ported in B09 using only the Cycle-8 Chandra data along with 
the XMM data, and thus the statistical significance of the line 
detection is now much larger. (We note that the measured line 
equivalent width of « 26 mA (Table |2]l is very close to, and 
consistent with, the value of 30 mA obtained by B09.) 

Following B09, we estimate the statistical significance of 
the WHIM and local lines using Monte Carlo simulations (see 
B09 for details). Briefly, we make two sets of simulated spec- 
tra for each observation. One includes both source and back- 
ground, and one is the background only. For the simulated 
background spectra, we used the model parameters obtained 
by fitting the real extracted background with a single pow- 
erlaw model. For the source spectrum, we use a model that 
includes the continuum and only one absorption feature, i.e.. 
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we include the Galactic line when estimating the significance 
for the WHIM line, and vice versa. We treat the two spectra 
in exactly the same way as the real data sets. In particular, 
as above we restricted the redshift of the WHIM line to lie 
between 0.028 — 0.032 to match the Sculptor Wall. 

First, we examined the significance of the Sculptor O VII 
line using only the Chandra data. In 100,000 simulations of 
the Chandra data without the Sculptor line, we obtained 64 
false detections with a reduction in the C-statistic at least as 
large as we measured from the real data; i.e., we conclude 
that the Sculptor O VII line is detected at a statistical signifi- 
cance of 99.936%, or 3.4(7, using only the Chandra data. This 
represents a substantial improvement over the 1.7ct detection 
we reported previously when using only the Cycle-8 Chandra 
data (B09). Second, when the XMM data are included, we 
obtain just 16 false detections out of 200,000 simulations cor- 
responding to a detection significance of 99.992%, or 4. Oct. 
Therefore, the Sculptor O VII Ka WHIM line is confirmed. 
If instead we use the standard data-weighted statistic in 
place of the C-statistic we obtain fully consistent results; i.e., 
we detect the WHIM line at a significance level of 3.3ct using 
only the Chandra data and at 4. Oct using both the Chandra and 
XMM data. We found these results are consistent with those 
obtained using C-statistic. 

The redshift of the WHIM line is well constrained sepa- 
rately for both the Chandra and XMM data (Table 2), and 
each are consistent with lying within the redshift range 0.028- 
0.032 representing the portion of the blazar's sight line in- 
tercepted by the Sculptor Wall. However, although the red- 
shifts are formally consistent at the « 90% confidence level, 
the Chandra data suggest a slightly lower value. (Note if the 
redshifts are tied between the Chandra and XMM data, then 
the WHIM line redshift is 0.0310!jJ:J]j;[!^ and the line is de- 
tected at a slightly lower significance level: we found 32 false- 
detections in 200,000 simulations, yielding a significance of 
99.984%, or 3.8ct.) A similar offset between Chandra and 
XMM is also observed for the local line. These differences 
reflect the approximately 40 mA wavelength shift between 
the XMM and Chandra data first reported in B09. The exact 
reason for the shift is unclear, since the typical XMM RGS 
wavelength accuracy is ^ 10 mA ; however, it is noted that 
occasionally the RGS wavelength scale can be in error by this 
magnitude possibly owing to operation issues (A. Rasmussen, 
private communication). 

For the column density and Zj-parameter we list their 90% 
errors on one interesting parameter in Table 3 and display 
contours representing their joint 68% and 90% errors on two 
interesting parameters in Figure |5] While the data do not 
place a strong constraint on b, we obtain a best-fitting value of 
96 km/s and a 90% upper limit of 455 km/s, which is an im- 
provement over our previous study in B09. The constraints on 
b are consistent with theoretical expectations for the WHIM; 
e.g., we expect even with shocks and turbulence the velocity 
should not be much higher than 500 - 600 km s"' in these 
superstructures (e.g., Kang et al. 2005). 

As shown in Figure |5] the constraint on the O VII col- 
umn density is sensitive to b. For > 150 km/s, the col- 
umn density lies in a fairly narrow range with a typical value, 
Novii ~ 2 X 10'^ cm"^, that is fully consistent with cosmolog- 
ical simulations. However, for smaller b the column density 
is not as well constrained: for b<5Q km/s the column density 
can take a large range of values from below Novii ~ 10'^ cm~^ 
to values exceeding Nqvii ~ 10'^ cm"^. Since O VII column 
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Fig. 5. — Displayed are the 68% and 90% confidence contours on two in- 
teresting parameters (AC = 2.30,4.61) for the Sculptor Wall absorption line. 
The (red) cross is the best fit while the (green) hatched area indicates the re- 
gion disfavored by cosmological simulations and observations of local X-ray 
absorption lines (see text). 

densities exceeding Novii ~ 10'^ cm"^ are neither easily pro- 
duced in cosmological simulations of the WHIM (e.g.. Fang, 
Canizares, & Bryan 2002) nor are they measured for the lo- 
cal X-ray absorption lines (e.g., Williams et al. 2005, 2007; 
Wang et al. 2005; Fang et al. 2006), we indicate such disfa- 
vored values of large column density with a (green) hatched 
region in Figure |5] 

Following the same procedure we also estimate the detec- 
tion significance of the Galactic O VII absorption line. This 
line was detected in the observations reported in B09, with a 
joint significance of about 3ct. With the new observation, the 
reduction of the C-statistic for the Chandra + XMM data is 
9.6. For the simultaneous fit of the Chandra and XMM data 
we obtained 145 false detections in 10,000 Monte Carlo tri- 
als for a 2.5ct detection significance. Using only the Chandra 
data, we obtained 1607 false detections for a 1.4ct detection 
significance. The somewhat smaller detection significance of 
the local line we have obtained with Chandra compared to 
B09 (who obtained 1 .7ct) can be understood in terms of the 
line equivalent width. Here we have measured an equivalent 
width of « 17 mA (Table |2]i for the local line, whereas B09 
obtained « 30 mA with a Ict error of 1 1 mA. While the val- 
ues are formally consistent at the 1 .2ct level, the best-fitting 
value has been cut almost in half, implying the line is weaker 
than previously estimated (note the O VII column density for 
the local line is very poorly constrained and, while the best- 
fitting value is larger than that of the Sculptor line, the lower 
limit is even lower.) 

4. SYSTEMATICS 
4.1. Anomalous Exposure #10498 
For one of the exposures (obs. #10498), there is a significant 
absorption feature located near 22.05 A which does not cor- 
respond to a known instrumental feature, nor is it present in 
any of the other Chandra exposures or the XMM observation. 
There is no known large-scale galaxy concentration along the 
blazar's sight-line at the redshift (z ^ 0.02) corresponding to 
the feature. One possibility is that the feature represents the 
z ^ 0.03 Sculptor line which has been shifted (by ^ 250 mA) 
as the result of a temporary gain non-linearity in the HRC-S. 
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Fig. 6. — Left panel: Background spectra for obs.#10577 between 20 and 25 A. Red shows result from the standard extraction procedure, and black is from 
using the improved pulse-height filter. Right panel: source spectrum of obs. #10577 between 15 and 25 A. Red line shows the model convolved by all the orders 
(first to the sixth orders), and green line shows only the first order contribution. 



We believe this is unlikely because a gain non-linearity should 
produce a shift no larger than 100 mA(Chung et al. 2004), 
and there is also no corresponding evidence for a shift in the 
local line. Instead, the transient feature is more likely intrin- 
sic to the blazar, a possibility we discuss in more detail in a 
separate paper (Fang et al. 2009, in preparation). 

Transient absorption line features were first discovered in 
early X-ray observations of blazars, and often interpreted as 
intrinsic to the source and can be characterized by, e.g., an 
absorption edge or a P-Cygni profile (e.g., Canizares & Kru- 
per 1984; Krolick et al. 1985; Madejski et al. 1991). These 
absorption profiles are sufficiently broad so as to impact sig- 
nificantly the local continuum shape near the Sculptor WHIM 
line, complicating interpretation of the WHIM line in this sin- 
gle exposure. For example, if the transient line shape follows 
a P-Cygni profile, it would contribute extra emission in the re- 
gion of the Sculptor WHIM absorption line which could then 
be compensated for by a stronger, more significant, WHIM 
line. We defer a detailed investigation of the modeling of the 
transient line to our dedicated paper on the transient feature. 
Here we estimate the maximum reduction in the WHIM sig- 
nificance we might expect when including exposure #10498. 

Such a conservative estimate is provided by assuming the 
transient line is Voigt-broadened like the WHIM line, al- 
though it is unlikely that a line intrinsic to the blazar will 
be thermally broadened. In this case, the transient absorp- 
tion line contaminates the 22.25 A WHIM line by contribut- 
ing excess absorption, lessening the need for the WHIM line. 
If we include exposure #10498 in our simultaneous Chandra 
and XMM fits while treating the feature as a separate, Voigt- 
broadened absorption line present only in this one data set, 
then we obtain fully consistent results for the WHIM absorp- 
tion line parameters, though the statistical significance of the 
line is slightly reduced (99.986%, or 3.8(t). Similarly, the sta- 
tistical significance of the WHIM line decreases from 3.4a 
to 3.2(7 (99.840%) when using only the Chandra data. Given 
such systematic uncertainties owing to the specific treatment 
of the transient line shape, and since this conservative esti- 
mate of the reduction in the WHIM significance does not have 
a qualitative impact on our finding of a clear WHIM detection. 



we exclude exposure #10498 from our analysis. 

4.2. Background 

Due to the malfunction of the anti-coincidence shield, 
the HRC-S has a higher background than expected prior to 
launch. The standard procedure to estimate the background 
involves defining a special region on the HRC-S to extract 
a background spectrum, and then perform filtering based on 
the photon energjQ Recently, the CXO has refined the 
background estimation procedure with an improved pulse- 
height filter that also includes time-dependent gain correc- 
tions (Wargelin et al. 2009). As in B09, we have adopted 
the improved procedure because, indeed, we find it delivers a 
superior background reduction. 

In the left panel of Figure|6]we show the background spec- 
tra for obs. #10577 between 20 and 25 A. Red shows the result 
from the standard extraction procedure, and black is from us- 
ing the improved pulse-height filter. The new filter reduces 
the background rate by almost 50%, and the variation is 
also significantly reduced. 

The smaller background produced by the improved proce- 
dure also leads to a higher detection significance for the O VII 
line. When preparing the Chandra spectra using the previous 
background estimation procedure without the improvements, 
we obtain a slightly lower detection significance (3.6a) when 
adding the Sculptor line. 

4.3. Higher Order Contribution 

As noted previously, the HRC-S cannot separate photons 
arising from different diffraction orders. The result is that 
the higher-order photons act as an additional "background" 
component on top of the dominant first-order spectrum. In 
this section we explore the sensitivity of our analysis to our 
model of the higher orders by taking the extreme alternative 
of ignoring them altogether. 

Taking obs.#10577 as an example, we estimate the contin- 
uum made by the first order photons. We convolve the model 
spectrum (obtained by fitting the source spectrum with the to- 
tal response matrix) with the first-order response. The result is 

see ,http://cxc.harvard.edu/cal/Letg/Hrc_bi7] 
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shown as the green curve in the right panel of Figure |6] Over 
21 - 22.5 A, the first-order contributes about 80 - 85% of the 
total photons. 

To assess the impact of ignoring the noise contributed by 
the higher-order photons in our analysis, we have performed 
our spectral fitting using the RMF and ARF files contain- 
ing only information on the first-order photons. Furthermore, 
since the first-order RMF does not account for the diffraction 
of photons with energies outside of the 21 - 22.5 A band of 
interest, we do not need to restrict the power-law index of 
the continuum models for each observation. In this way we 
are able to assess how much the detection significance of the 
O VII line is degraded using a higher background and a more 
flexible continuum model. Despite these changes, we find that 
the C-statistic is still reduced by 19.8, very close to the value 
we obtained when properly accounting for the higher-order 
photons. 

We also note a weak correlation (not shown) between the 
slope of the power-law index for the broad-band continuum 
and the reduction in the C-statistic obtained for each Chandra 
exposure when including the Sculptor line. This behavior is 
expected because a flatter continuum implies a larger propor- 
tion of higher-order photons from shorter wavelengths being 
diffracted into the 21 - 22.5 A band, thus effectively lower- 
ing the S/N. The trend is fairly weak because the exposure 
time and flux of each observation also affect the S/N, and the 
variation in photon index between the observations is modest. 

4.4. Choice of Spectral Binning 

In our analysis we have adopted a conventional approach 
where the individual spectra have been rebinned so that each 
bin has counts above some minimum threshold value. The 
thresholds adopted (40 for Chandra, 75 for XMM) were cho- 
sen to be consistent with B09. However, other binning choices 
are possible which can yield different results for the statistical 
significance of the WHIM line (e.g., Gumble 1943, Humphrey 
et al. 2009). To assess the magnitude of possible variations 
in the WHIM line significance arising from the binning pro- 
cedure, we considered a more complicated scheme, focusing 
only on the Chandra data for simplicity. 

In this scheme each Chandra spectrum was first analyzed 
separately for a suite of different binning definitions (consid- 
ering different rebinning factors and shifts in the location of 
the first bin). For each individual spectrum we selected the 
binning definition that gave the maximum C-statistic differ- 
ence obtained from fitting models with and without the Sculp- 
tor WHIM line only to that individual spectrum. After the 
binning definitions were determined in this way, we computed 
the WHIM line significance by first simultaneously fitting all 
of the Chandra spectra to determine the total change in the 
C-statistic. This change in the C-statistic was then compared 
to those obtained from Monte Carlo simulations as before. In 
particular, as above we restricted the redshift of the WHIM 
fine to lie between 0.028 — 0.032 to match the Sculptor Wall. 
In this case, however, the more complicated rebinning pro- 
cedure is also incorporated in the Monte Carlo simulations; 
i.e., the binning of each mock spectrum is determined spe- 
cially for each Monte Carlo simulation as specified above. 
We find that when enacting this binning procedure, we ob- 
tain a statistical significance of 3.2<t for the WHIM line (132 
false detections over 100,000 simulations), similar to that ob- 
tained (3.4(t) with the conventional binning approach. (The 
quoted significance estimates are accurate to w O.lu consid- 



ering Poisson noise in the number of simulations/trials.) 

5. DISCUSSION & SUMMARY 

Using a deep, 500 ks follow-up Chandra observation, we 
confirm at higher significance the O VII Ka absorption line 
in the Sculptor Wall (z ^ 0.03) in the sight-line of the blazar 
H 2356-309 (z = 0.165) first reported in B09. We detect the 
O VII line at the 3.4a significance level when using only the 
Chandra data, a substantial improvement over the 1.7a level 
reported previously when using only the Chandra Cycle-8 
data (B09). The significance increases to 4. Oct when includ- 
ing the existing XMM RGS observation in a joint fit with the 
Chandra exposures. While the significance of the Sculptor 
line improves substantially with the addition of the new Chan- 
dra data, the significance of the local (z ^ 0) O VII line is 
slightly reduced compared to that found by B09 {2.5a vs. 3a). 
Finally, in one of the Chandra exposures we find evidence for 
a transient absorption line intrinsic to the blazar 04. lb . which 
we discuss further in a separate paper (Fang et al. 2010, in 
preparation). 

5.1. Interpreting the Sculptor line as Intrinsic to the Blazar 

While the possibility that the WHIM line is instead intrin- 
sic to the blazar cannot be categorically ruled out with the 
current data, we believe this to be highly unlikely. The red- 
shift of H 2356-309, determined by stellar absorption fea- 
tures from the host galaxy is z = 0. 165 ± 0.002 (Falomo 199 1 ). 
Therefore, the most likely candidates for an intrinsic absorp- 
tion line would be O VII Ka with an implied outflow veloc- 
ity of w 0.12c, or O VIII Ly-a, inflowing at -2000 km 
(in contrast to the expected outflows in the polar region; e.g. 
Proga 2005). In the latter case, if we consider that the line 
arises from a cloud of inflowing material, and if we make the 
reasonable assumption that the inflow speed does not exceed 
the freefall velocity, it follows that the cloud must lie within 
D — 6.4 X lO'^Ms cm of the black hole (which has mass 
Mg X IO^Mq). As there are no significant absorption lines 
of other species at the same redshift, we assume that O VIII is 
close to its peak ionization fraction (— 0.5), which occurs for 
ionization parameters, ^ = L/D^nu w20^0 erg cm s"' (Kail- 
man & Bautista 2001), where L (erg s"') is the source lumi- 
nosity and nu is the number density of hydrogen in the cloud. 
Given the luminosity of the AGN (L — 5 x 10"*^ erg s"', as- 
suming isotropic emission), maintaining this ^ would require, 
hh > 3 X lO'^M^^ cm"-'. Assuming a uniform, spherical cloud 
with Solar abundances, the measured O VIII column density 
(— 10'^ cm"^) would thus require a size, /?cioud = A^h/wh < 
10'^M| cm. The persistence of the WHIM absorption line 
over the time f 1 .5 yr spanned by the XMM and Chandra 
observations would therefore imply a velocity perpendicular 
to the line of sight, Rdoui/t < 0.2Mg km s"'. Even allow- 
ing for black holes as massive as IO'^Mq, this would require 
an extremely fortuitous (and improbable) orientation of the 
cloud's velocity vector to the line of sight. 

Interpreting the line as outflowing material is also problem- 
atical. While high-velocity outflows have been reported in 
some other AGNs (see, e.g.. Pounds & Reeves 2009), typi- 
cally they exhibit features produced by several different ion 
species. Although we cannot rule out there being other, weak 
lines, they are certainly not as prominent as the implied O VII 
line. This situation is to be contrasted with the interpreta- 
tion as a WHIM line, where it is Ukely a coUisionally ionized 
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plasma in which a dominant O VII Ka line is expected. Fur- 
ther observations of H 2356-309 will be useful to investigate 
whether weaker lines are present. Nevertheless, reconciling 
such a putative outflow with the appearance of a transient line 
at ^22.05A (Section l^Tt would require the outflowing mate- 
rial to exhibit a wide range of line-of-sight velocities. 

5.2. Density of the Sculptor WHIM 

Our clear detection of the WHIM absorption line allows us 
to make a rough estimate of the three-dimensional density of 
the WHIM gas. We obtain a 90% lower limit on the O VII col- 
umn density of 0.8 x 10'* cm"^ and a 90% upper limit on the 
Doppler b parameter of 455 km s"', each of which are con- 
sistent with cosmological simulations of the WHIM 03.2b . 
Assuming the absorber is uniformly distributed along the por- 
tion of the blazar's sight-line that intercepts the Sculptor Wall 
(- 15 Mpc from z = 0.028 -0.032)[3, the density of hydrogen 
is, 



„ = 6 X 10"''(cm-^) 
Nov II 



2 X 101'' cm-2 I 1 



fo 



Q.lZo 



15 Mpc/ 



where we have used a typical value of Novii corresponding 
to b > 150 km where the inferred column density is only 
weakly dependent on b (Figure |5]l, fovii is the O VII ioniza- 
tion fraction, L is the path length, and we assume a solar ratio 
of O/H and a metallicity, Z, of 0. 1 solar typical of low-density 
gas produced in cosmological simulations (Dave et al. 2001; 
Cen & Ostriker 2006). We assume fovii ^ 1 because the data 
do not show evidence of a very significant, redshifted O VIII 
Ka line at --19.55 A. (The study of the lines other than O VII 
Ka will be the subject of a future paper) Hence, this density 
in equation (1) translates to a baryon over-density of 6 ^ 30, 
which is consistent with the peak of the WHIM mass frac- 
tion predicted by cosmological simulations (Dave et al. 2001; 
Cen & Ostriker 2006). Moreover, the cosmological simula- 
tions predict that WHIM gas at this over-density should have 
a temperature of lO*" K. This temperature also happens to 
correspond to the peak in the O Vll ionization fraction, con- 
sistent with our assumptions. 

Since the Sculptor line provides key new evidence for 
WHIM with over-density and temperature typical of (> 50%) 
of the WHIM gas, it also provides a crucial link between 
the results obtained from O VI absorption and X-ray emis- 
sion studies of the WHIM. First, the O VI absorbers typically 
probe the low-density, low-temperature portion of the WHIM, 
because the O VI ionization fraction peaks near T ^ 3 x 10^ 
K assuming collisional ionization. If the O VI absorbers are 
produced by photo-ionized gas, then the gas is even colder 
(T < 10^ K). In all, probably 10-20% of the WHIM gas is 
probed by the O VI absorbers (Danforth & Shull 2007; Tripp 
et al. 2008; Thom & Chen 2008). On the other hand, the 
WHIM gas probed by X-ray emission studies typically has a 
higher temperature (~ 10^ K) and higher over-density 150; 
e.g., Werner et al. 2008), which represents the higher-density 
tail of the WHIM distribution (Dave et al. 2001; Cen & Os- 
ti-iker 2006). 

We note that while the Hubble flow could contribute to the line broad- 
ening, the effect is very uncertain in a superstructure like the Sculptor Wall, 
where local gravitational interactions between galaxies may dominate the 
global Hubble flow. 



5.3. Future Work 

Due to the Umitations in the spectral resolution of the Chan- 
dra and XMM data, we are unable to map in detail the spatial 
distribution of the absorber along the line-of-sight. Because 
the Sculptor Wall is densely populated with galaxies, it is pos- 
sible that the O VII Ka WHIM absorber we have detected 
actually corresponds to the halo of an individual galaxy or 
group, rather than to the entire ^^15 Mpc sight line. Indeed, 
such X-ray absorption is seen in our Chandra and XMM spec- 
tra (i.e., the z = O VII line) and has been detected repeatedly 
in and/or around the Milky Way (e.g., Nicastro et al. 2002; 
Fang et al. 2003; Rasmussen et al. 2003; Yao & Wang 2004; 
Williams et al. 2005). Although it is still under debate whether 
these absorption lines are produced by hot gas in the disk, in 
the halo of the Milky Way, in the intra-group medium of the 
Local Group, or some combination of the three, several lines 
of evidence (e.g., Wang et al. 2005; Fang et al. 2006) do sug- 
gest the hot gas in our Galaxy may be responsible for absorp- 
tion. 

One method to address the spatial distribution of the Sculp- 
tor WHIM line is through X-ray imaging of the WHIM emis- 
sion. If the Sculptor WHIM actually arises from the halo 
of a MiUcy-Way-sized galaxy, there should be no significant 
WHIM emission outside the virial radius, corresponding to 
about 6' at z ~ 0.03. Although imaging studies of the WHIM 
are better suited to much higher over-densities than we infer 
for the Sculptor Wall gas, it should nevertheless be possible 
to detect it in emission using offset pointings at least 6' away 
from the blazar with Suzaku and XMM-Newton, if it is actu- 
ally extended throughout the entire Sculptor Wall. 

Our successful detection of the O VII Sculptor WHIM ab- 
sorption line demonstrates the viability of our observation 
strategy using Chandra and XMM; i.e., deep (0.5-1 Ms) 
targeted observations of a modest number (5-10) of bright 
blazars behind known large-scale galaxy structures is an ef- 
fective means to detect and study the WHIM in X-ray ab- 
sorption with current observatories and is much more effi- 
cient than blind searches. To achieve a complete census of the 
WHIM gas in the X-ray band, higher quality data are needed 
such as would be provided by the planned International X-ray 
Observatory (IXOH The grating for the IXO has an effective 
area of '--^ 1 , 000 cm"^ at the soft energy band, two orders-of- 
magnitude larger than that of Chandra gratings. This means 
a few ten ks observations of H 2356-309 with the IXO would 
achieve the same level of detection as we reported in this pa- 
per Also, the much higher spectral resolution of the IXO grat- 
ing will clearly resolve the WHIM line structure, therefore 
providing rich details of the WHIM gas. Finally, such tele- 
scopes will make the "blind"-search strategy possible since it 
can probe the WHIM gas in very low density filaments, and 
complete the census of the WHIM gas. 
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